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THE PASSAGE OP AN INFINITE SWEPT AIRFOIL 


THROUGH AN OBLIQUE GUST 

By John J . Adame zyk 
United Aircraft Research lahoratories 


SUMMARY 


An analysis is presented which yields an approximate solution for the 
unsteady aerodynamic response of an infinite swept wing encountering a . ' 

vertical oblique gust in a compressible stream. The approximate expressions 
are of closed foirni and do not require excessive computer storage or computation 
time, and further, they are in good agreement with the results of exact theory. 
This analysis is \ised to predict the unsteady aerodynamic response of a 
helicopter rotor blade encountering the trailing vortex from a previous 
blade. Significant effects of three-dimensionality and compressibility are 
evident in the results obtained . In addition, an approximate solution for 
the unsteady aerodynamic forces associated with the pitching or plunging 
motion of a two-dimensional airfoil in a subsonic stream is presented. The 
mathematical form of this solution approaches the incompressible solution as 
the Mach number vanishes, the linear transonic solution as the Mach number 
approaches one, and the solution predicted by piston theory as the reduced 
frequency becomes large . 



INTRODUCTION 


The prediction of the unsteady aerodynamic response of an airfoil, to a,, 
vertical gust velocity field has long "been of interest to aeroelast.icians and 
acousticians. For the most part aeroelasticians have used the incompressible 
two-dimensional theories of Kussner (Ref. l) and Sears (Ref. 2) to predict 
the unsteady response. In estimating the aco-ustic field generated by airfoil- 
gust interactions the acoustician typically, uses the unsteady lift obtained 
from one of the previous theories to determine the strength of the acoustic 
dipole source which is then assmed to replace the airfoil in the flow field 
(e.g.. Curie in Ref. 3 )- 

lypical of the many unsteady aerodynamic problems in which these 
incompressible theories have been used in the past are: (l) rotor blades 

passing through the wakes of stator blades in turbomachinery, (2) an airfoil 
interacting with a turbulent gust, and (3) helicopter rotor blades encounter- 
ing the tip vortices from preceding blades . In each of these problems some 
uncertainty arises as to the applicability of two-dimensional incompressible 
aerodynamics. It is not surprising, therefore, that several analytical 
studies have appeared recently which treat the complexities of three- 
dimensionality and compressibility in the unsteady problem. For example, 
Filotas (Ref. k) obtained a closed form solution for an oblique sinusoidal 
gust encountering an infinite wingj however, this work is limited to a 
two-dimensional airfoil of zero sweep (i.e., the incoming flow is normal 
to the leading edge line of the Wing) in an incompressible stream. Graham- 
(Ref. 5) included the effects of compressibility and three-dimensionality 
but neglected the effects of sweep. His analysis is based on a numerical 
solution of the governing differential equations. Adamczyk (Refs. 6 and T) 
also included the effects of compressibility and three-dimensionality, but 
the form of the solution was expressed in terms of an infinite series of 
Mathieu functions which are difficult to eval\iate analytically. Another 
recent analysis by Johnson (Ref. 8) included the effects of compressibility 
and three-dimensionality; however, his approach is tailored towards analyzing 
the response of an airfoil to a free rectilinear vortex. Hence, the 
researcher has had no simple compressible three-dimensional analogue of the 
Sears solution available to him and was required to resort to the more 
restricted theories of Filotas or Sears for predicting the aerodynamic 
response of an airfoil to a vertical gust. 


To overcome the deficiencies described above, an analysis has been 
developed which generates simple approximate expressions for the unsteady 
compressible lift and moment response functions for an infinite swept airfoil 
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encountering a three-dimensional oblique gust (see Fig. 1 and Fig. 2). The 
approximate expressions are of closed form and do not require excessive 
computer storage or computational time . This report describes the development 
of these approximate expressions, and presents a comparison of results with 
an exact theory and sample computations of the unsteady loading of a 
helicopter blade encountering a vortex. 
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NOMENCLATURE 


a 

A 

n 
b , 

B 

n 

C 

ce 

n 



conj 

E 

f 


gj_, gg 
h 


H H 


H 


^0' 

^ 0 ’ ■^1 


n 


K 


/ 

speed of sound 

. T 

coefficients in the series for defined in Eq. (46) 

semichord of the airfoil ' - 

( 2 ) 

coefficients in the series for p 

cosine Fresnel Integral 

symmetric Mathieu functions 

lift and moment coefficients per unit span 

complex conjiagate operator 

complex Fresnel Integral, defined by Eq^ ( 85 ) • ■ 

function defined by Eq. ( 96 ) or Eq. ( 151 ) 

functions defined by Eqs. (l35)^ (l36)p'(l38), and (l39) 

displacement of the airfoil from the x^ plane, displacement 
of the vortex from the x^ plane 

amplitude of plunging motion .. . 

Hankel functions, first kind, orders 0 and 1 

complex constant, defined by Eq. (31 modified acceleration 
potential, defined by Eq. (93) 

integral defined by Eq. (99) 

modified Bessel functions, first kind, order 0 and 1 
Bessel functions, first kind> orders 0 and 1 
Bessel •.'’unction> first kind, order h 
reduced frequency 
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K 

gust wave number 

m 

integer 

M 

Jfech number of the free stream, U/a, function defined by 
Eq. (97) 

M* 

Mach number of the gust relative to the oncoming flow, 
defined by Eq. (l?) 

Me 

n 

Mathieu-Hankel function, second kind, associated with ce 

' n 

n 

(3 M (^ ) ^ £ ^ 

n 

integer 


incompressible pressure distrihution, defined by Eq. (TO) 

P > P 

s 2 

first and second terms of the solution for p for large 7 

P 

perturbation pressure 

Pq ^ ^ 

: steady pressure, associated with, the flow U 

- 

:= 

r 

radial distance from the midchord line of the airfoil, 
distance measured in the plane of X , X , defined by 
Eq. (i4t) ^ 

R 

dimensionless radial distance 

S 

Sears function, or sine Fresnel Integral 

sgn 

sign operator 

t . ■. 

time 

f 

nondlmensional time, defined by Eq. ,(43) 

\ 

lift transfer function, defined in Eq.' .(22) 

\ 

moment transfer function, defined in Eq. (23) 
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integration variable 

free stream velocity j 

component of perturbation velocity 
complex amplitude of Ug 
Cartesian coordinates, Pig. 1 

dimensionless coordinates, defined by Eqs. (ll), (12), and (13 
variable 

complex constant, defined by Eq. (25) 
angle of encounter. Fig. (l) 
constant, defined by Eq. (43) 
amplitude of pitching motion 
frequency parameter, defined by, Eq. (l4) 
cutoff parameter, defined by Eq. (l?) 
circulation of free vortex 
constant defined by Eq. (79) 
functions defined by Eqs. (13T) and (l40) 
angle of sweep. Fig. 1 

elliptic-hyperbolic coordinates, defined by Eqs. (46) and (47) 
density fluctuations 
mean density 
integration variable 

perturbation velocity potential defined in Eq. (6) 
complex amplitude of 0 



complex amplitude of 4> for an incompressible fluid 


'’’iNC 

4 ,( 1 ) 

4 ,( 2 ) 


noncirculatory component of $ 
circulatory component of 4> 
Laplaclan operator 


Superscripts. 


I 

S 


A 


incident disturbance 
scattered field 

Fourier transform of a variable 


Subscript 

1/4 


quarter- chord 



FORMULATION OF THE MATHEMATICAL FROBLEM 


Governing Eq\aation 


Figiore 1 shows the geometry of the problem. An infinite swept airfoil 

of chord 2b is held rigidly in a compressible fluid which is flowing over 

it with uniform velocities Ucos0 in the direction, and Usin^ in the 

direction, Superposed on U (see Fig. 2) is an unsteady three-dimensional 

disturbance Uj[. The pressure and density are pq + p, and Pq +P, respectively> 

where the subscripted part denotes the^mean value, while the remaining part 

is associated with the disturbance. Far upstream of the airfoil, the 

ITT 

perturbations Uj;, p, and p are given functions, denoted by Uj^ , p , p . 

The problem is described by the customary fluid dynamic equations for 
a compressible flow (see Ref. 9)* The equations of conservation of mass 
and momentum are linearized by assuming that the perturbations are small 
compared to the mean values thereby allowing the squares and products of 
the perturbations (and their derivatives ) to be neglected while first-order 
terms are retained. Hence, 


dp ^ dp . ^ dp duj 

— ^ -I- U cos 9 + U sing- + T — - -0 


( 1 ) 


dt 


dx 


dXj dXj 


I doi ^ du\ . - \ 


dp 

^ -d^ 


( 2 ) 


For small fluctuations, the change in pressure is proportional to the change 
in density. 


p = a^p 


(3) 


where a is the speed of sound. 
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Far from the airfoil the pressure, density, and velocity field must 
eq\jal the incident disturbance. This condition 




( 4 ) 


is most conveniently satisfied by taking 


{piyo .Mi} = 


(5) 


and imposing the condition that p®, p® and u^^® vanish far from the airfoil. 

The quantities p^, p^ and u^^ must satisfy Eqs . (l), (2) and (3) because 
far from the airfoil they constitute the entire disturbance. Hence, because 
the equations are linear, the quantities p®, p®, and Uj[® must also satisfy 
Eqs. (l), (2), and (3). Furthermore, since the fluid is inviscid, no change 
in vorticity can take place as the fluid passes over the. airfoil. Hence, Uj^® 
niust be an irrotational field, and may be represented by the derivative of a 
velocity potential 



d<f> 

dxi 


( 6 ) 


The momentum equation (2) provides a relation between the pressure and the 
potential, namely, ’ • 


r^S = 


/ Atb 


(T) 


and when this and the equation of state (3) are substituted into the 
continuity equation (l), one obtains the governing differential equation 
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for the problem 


(■ 


2 2 2 2 d if> 

M COS 0 r -M sin 9 


dxt^ 


‘I 


ax,2 


) 


2 /Mcos^ a> Ms\n9 

* Q dx,at ''' a dxjdt 


+ M COS 0 sin 0 


d^4> 




d% 

— ? = 0 

af2 


( 8 ) 


Here M = U/a is the Mach numher of the mean flow. 


Boundary Conditions 

A solution of Eq. (8) is sotight which describes outgoing waves which 
must decay as T~i for large distance, r, from the airfoil (cf . Ref. 7 )- On 
the surface of the airfoil the velocity normal to this surface miist vanish. 
Hence, 


d<f> I 


axg 


U 2 X2 = 0 - b < X, < b 


( 9 ) 


U2^ being a given function of space and time. 

The present analysis assumes that the vertical velocity field induced 
by the incident disturbance on the airfoil is that of an oblique sinusoidal 
gust convected at the free stream velocity U. The mathematical form for 
this velocity field is 


^ 2 ' Ug 


|i[K cos aX| + Ksin a X 3 -K ucos ( 6 -a)t]| 


( 10 ) 


where the sweep angle 6 and the encounter angle a are defined in Figs . 1 and 
2. Since the gust is assumed to be convected with the mean flow velocity U, 
the pressure and density fluctuations associated with this velocity field 
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are zero (l.e«> = O). Dimensionless coordinates and a frequency 

parameter are Introduced for convenience in sulDsequent operations. ■ 

x,= x,/b ( 11 ) 


I 

X2 = X2(i-m^cos^ 5)^ /b 


( 12 ) 


Xg= Xj{ I - M^COS^ 5)^ / b 


( 13 ) 


“R b cos g 
l-M*COS*5 


( 1 ^) 


The solution for the velocity field Uj^® must have the same frequency and 
the same spanwlse wave number as the incident gust. In addition it is found 
convenient to Include an exponential factor in the X]_ direction of the potentla, 
Thus the following form for <f> is chosen: 




4 > = «>(X|,X2)exP 


^cos^5x, - Ku cos(6-a)t + 


Kb sin g 



( 15 ) 
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When the substitutions (ll) through (15) are made in the governing eq\aation, 
( 8 ), one finds that 4 > must satisfy 


\ 



d^<I> 2 

5- 1 y ♦ = 0 

dX2 


\ 

where y is a compressible cutoff parameter defined as 


K b sin g 



K b sin aV^*z I 
I -M^cos^ff 


(16) 



(IT) 


Note that the sign of y^ is dependent on whether the parameter M* 
is greater than or less than one . This parameter is proportional to the phase 
velocity of the disturbance along the span of the airfoil relative to the 
oncoming flow (see Fig. 3 )* When M* is greater than one, the relative phase > 
velocity of the disturbance is supersonic. However, if M* is less than one, 
the relative phase velocity is subsonic. The behavior of the solution of 
Eq. (16) at large distances from the airfoil is critically dependent on the 
value of M*. This dependence can be shown by examining the asymptotic limit 
of the solution of Eq. (16) for an isolated acoustical source (cf. Ref. 10 ). 
This solution has the form 


/-r — r 

* ~ R=VX, +X2‘" -‘00 (18) 


Eq-uation (18) requires <J) to decay exponentially with distance for imaginary 
y which, from Eq. ( 1 T), Is equivalent to M*< 1 . If the cutoff parameter y 
is real (l.e., M*> l), Eq. (18) is the asymptotic form for a cylindrical 
acoustical wave propagating outward from the origin at R •» 0 . 

The boundary conditions on are 


d<P 

3X2 




U2e'^’‘‘ for -i< x,< l,X2=0 


(19: 


12 



which was derived from Eq. (9) and 


_ 

\ <J> 2 Q5 R ^ 00 

» ' . 


( 20 ) 


which is the asymptotic radiation condition of Eq. (18). In addition to the 
hovmdary conditions in Eq. (19) and (20) a Kutta-Jovikowski condition must he 
imposed at the trailing edge of the airfoil. This condition requires that 
the presstire Jump across the airfoil vanish at the trailing edge. This 
requirement can he mathematically expressed in terms of 4> and its derivatives 
hy substituting Eq. (15) into Eq. (7 ),to yield 




U COS© 

b 


axi 



exp 


|-i/3M^cos^©X| 


^ iKb.sinaX3_ 
y I - M^cos^© 


( 21 ) 


which equals 


p®= 0 at X, = I, Xg = 0 


( 22 ) 


The governing differential equation' (16) and the boundary conditions, 

Eq. (19), (20) and (22), form a boundary value problem for $ whose solution 
is dependent on only two variables, the cutoff parameter y, and the frequency 
parameter y©. The construction of the approximate solutions to this bou. dary 
value problem was divided into three parts : (l) a solution valid for small 

positive values of 7, (2) a solution valid for negative values of 7^, and (3) 
a solution valid for large values of 7^, each of which will be discussed beibw 
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SOLUTION OP THE BOUNDARY PROBLEM 


General Considerations 



An important objective of this analysis is to determine the .unsteady I 
lift and moment coefficients (per unit span) resulting from the. encounter: of • 
an infinite swept airfoil with a three-dimensional oblique sinusoidal gust! inr 
a compressible stream. This is accomplished by integrating the zero and first 
order moments of the pressure distribution over the chord of the airfoil; 
however, the pressure distribution is a highly specialized function of the 
cutoff parameter, and, as shown below, the solution changes its character 
for various regimes of 7. In general, the lift and moment' coefficients can 
be written in the form 


! 27rlj2Cos 9 

u -/ 1 - M^COS^ 


Tl6 


i K sin aX3~UK cos(0-Q)t] 


(23) 


'M - 


7T Tig COS 6 

2U>/|-m2cos^0 


i[R sina X3 - UK cos(0-o)t] 

Uw e , 


(2U) 


where! Tl and Tj^j are the lift and moment transfer, functions, and the expressions 
within the brackets | | represent the quasi-steady lift and moment .s response . 
The functions T^ and TJ^,J will now be derived for. the various values of 7^ of 
interest here. ' 


Solution for Small Values of 72:0:. 

The approximate solution of Eq. ( 16 ) for small positive values of 7 was 
obtained by expanding the exact solution to this problem, taken from Ref. J, 
in a power. series of y. The -analysis- presented in Ref. V^was developed , by. 
linearly separating 4 > into two components. The first component, 
accounts for the noncirculatory flow field surrounding the airfoil. This 
solution for satisfies Eq. ( 16 ) along with boundary conditions in 

Eqs . ( 19 ) and (20) and has the following properties : (i) the circulation 

around the airfoil is zero, and (ii) singularities in pressure and velocity 
occur at the leading and, trailing edges , of the airfoil. The desired solution 
for is obtained by adding to O' a second component, 0'^', which is the 


ll 


circulatory solution. This new term cancels the trailing edge singularity 
of $vl; (i.e., causes the solution f or to satisfy the Kutta-Joukowski 

condition, Eq. (22)), and has its normal gradient vanishing on the airfoil. 

The details of the approximate solution for small positive values of 7 
will' he • found in Appendix I in which the final form for the pressure distri- 
bution on the surface of the airfoil is shown to he : 


s A>qUU;Cos8 /TTT 


I 


-iZJ|(/ 9 ) 

+ z 


exp 


-/Sm^cos^®X|+ ■ X3~^dcos{d-a)t| 


(I- 


(25) 


- I <X, < I 


where the variable 


iHo\/3) + i(y//3)^|Ho^'\y/2)-Ho^'\/3) 


h‘”(/8) 


(26) 


The equation for the lift and moment transfer functions Tj^ and Tj^, respectively, 
can he obtained by integrating the zero and first order moments of the pressure 
distribution over the chord of the airfoil. The resulting expression for the 
lift transfer function Tj^ is : 


Tl 


. Jo(^-iZJ|() 8 ) 


l+Z 


[jo(/3M^cos^e)+ iJ|(/8M^cos^0)]e 


ibKcosa 


(2T) 


while that for the moment transfer function T2^{ about the qiiarter-chord is : 

__ r - / . .1 

ibKCOS 


Tm ,/4 = e 


>sa ~ i Z J| [0) 2 2^\ 

\ + 2 -J(/^M^COS^0) 


( 28 ) 


+ ^^{0 ) - i 2 J| i^M^COS^d) + Tl 


2 2 , 
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Equations (25), (27), and (28) are used to detennlne the aerodynamic loading 
for small positive values of y. 


Limiting Behavior as y-^O'' 


■ I 5 : 


The influence of the cutoff parameter 7 on the lift transfer function ’ 
for small values of this parameter may he' shown ‘ty examining the limiting* 
form of Eq. (27) as y-^O"^. This limit is given by the equation ' = 


Tl = 


conj)s(/?)} - Jo^/?M^cos^8) 


• ij ^ 

xy T.- > 3 

- J . <i I'j 


+ iJ,(ySM^COS^ 


a)]e‘ 


iKbcos a 




(29) 


where conj | j denotes the complex conjugate of the enclosed qviantity, and 
where 


y*,„ j; 

^'"4)8 


b K [( M COS sin^a] . 

vfMcos8)^-sin^a 

,09 In 

(l-M‘^cos e) cos a 

4 cos a 


(30) 


can be obtained by manipulating Eq. (l7)* Note that S()8) is the two-dimensional 
gust response function derived by Sears (Ref. 2). The influence of y is. shewn 
by Eq. (30) to be of first order in K and (Mcos ^ - sin^a, and its effect on 

the mcjment transfer function is of secondary importance compared with the; 
parameter ^V^cos^B. ' Ah expansion of Eq. (28) for small;yalues of j8^^,cos.^6, . 
yields the result that the moment transfer function is a llhear function of 
this parameter. This result implies that the moment transfer function about 
the quarter chord vanishes as (l) the gust wave number if approaches its 
steady-state limit of zero, or (2) the speed of sound of the fluid medium 
approaches infinity. 

Solution for Values of y^sd 

A power series expansion of the result presented in Ref. 7 for the 
circulatory pressure distribution on the surface of the airfoil for negative 
values of 7 ^ was derived. It is shown above- that this approach leads to 
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v 4 ry simple expressions for the aerodynamic response functions (i.e., Eqs.' 

(25), (27), and (28)) for small positive values of T. However, when the 

o 

expansion technique was applied to the problem for y it yielded solu- 

tions which contained definite integrals whose values had to he determined 
numerically. This form of a solution is not satisfactory hecaxise it does 
not^.lead tp, a .simplified analytical proced\ore. Hence, an alternative 
approach, was dpyeioped .based on Graham's similarity rules (Hef. 5) and the 
analysis developed by .Filotas (Ref. V) for an oblique sinusoidal gust 
encountering an infinite airfoil in an incompressible stream. This 
alternative approach led to an approximate expression for the pressiire 
distribution on the surface of the airfoil which can be applied over the 
entire range - : \ J > ^ - 

i ". ' 2 

The details of the approximate solution for negative values of 7 will 

be found in Appendix II in which the final fom of the pressvire distribution 
is shown to be 


_ /’oUUgCos 6 H{y0) /PxT e 

Vi-M^cos^ V i+x; 


-l<X<l 



,8 = tQn’V//3 (33) 
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The lift transfer function Tj^ is computed hy integrating the pressure 
distribution, Eq. (31 over the chord of the airfoil, yielding 




( 3 ^^) 


• J 


The corresponding expression for the moment transfer function Tjj,j about the 
quarter-chord is 


Tm 


2 e 


iKb cos a 


In(r)4I,(X) 


Io(X+ iySM^COS^^) 


(35) 


+ I,(r+i/3M%s^0) 




r+i)3M^cos^0, 


H(7,J3) + Tl 


Equations ( 31 ), (3^), and (35) can be used. to compute the aerodynamic response 
functions over the entire range of negative V . 


Solution for Large Values of 1 7^1 

The results predicted by the analysis of Ref- 7 for the pressure distri- 
bution on the surface of an airfoil encountering a two-dimensional sinusoidal 
gust indicate that for large values of |7^| the pressure distribution appears 
to be composed of acoustical waves whose points of origin are the leading and 
trailing edges of the airfoil. The spatial influence of the acoustical wave 
emitted from the leading edge of the airfoil Increases as 7 increases, while 
the influence of the wave emitted from the trailing edge decreases . This is 
indicated by the behavior of the pressure phase angle (both theoretical and 
experimental) in Fig. 11 of Ref. 11 and is further corroborated by computa- 
tions performed during the development of the basic theoiy (Ref. 6). This 
behavior implies that for large values of ( 7^| a good approximation to the 
pressure field in the neighborhood of the leading edge (i.e., far upstream 
of the trailing edge) can be obtained by considering the • interaction of a 
g-ust with a semi^inf inite plate. For values of 7^ ^ 0 the pressure distri- 
bution p^^^ on the surface of the plate possesses a leading edge singularity, ' 
and decays with distance aft of the leading edge. Thus the pressiire distri-' • 
bution is a monotonically decreasing function which vanishes only at an.^ 
infinite distance downstream of the leading edge, and therefore the pressure 
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does not satisfy the Kutta-Joukowskl condition. This error 
ed hy adding to a pressure field p^^' which is eqml and 

opposite to p'^' downstream of the airfoil (i.e., 1, Xjg = O) and has' its 

normal gradient vanishing everywhere on the airfoil. 


distribution p 
can he correcte 


^The details of the approximate solution for large values of will he 
found in Appendix III in which the final form for the pressure distrihution 
on the svirface of the airfoil is shown to he 




(36) 


where 


^( 1 ) = 


/OqU02COS 6 


(i + i) 


*exp 


[yi-M^cos^® \/z)^fTr{0+r) 

{y- $m^cos^6)x( t iRsin a Xj- iuR cos (0-a)t| 


(37) 


and 


.( 2 )= 


1 Po ® 

-“-E^y2r(2-x;)| 



'-jSM^COS^0)x[-* 

■ iKsinax 3 -iUKcos (0-a)t 


exp 


(38) 


X,'=X,+I, -I < X|<1 

The mathematical form of the two-dimensional limit of Eq. (37) is that of an 
outward propagating acoustical wave whose point of origin is the leading 
edge, while -the equivalent limit of Eq. (38) yields a form which represents 
an outward' propagating wave whose apparent origin is the trailing edge . Once 
again the lift and moment transfer function are obtained by integrating the 
zero and first order moments of the pressure distrihution, Eq. (36), over. 
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the chord of the airfoil. The result for the lift transfer function Tj^' ■ , j 
■ is , 


T, -- 


I 


y->r(/3+y) 


(i + i) 


7T 


y-$ M^cos^® 


E \^z (y-/3M^cos^ej + [ i-e 


2i(y-/Si/co^0) 


2(y-j8M^cos^9) (y-/3M^cos^e) /-/SM^cos^S 


iE[y^] i yi7 

^ . ■ A .N 


2 2, 


g2i(y-/3M^COs?0) 

^y+/3M^cos^0 


2(/+ )9 M^cos^0 


while that for the moment transfer function about the quarter-chord is 

I /-Z Z-T-1 

E|yy.(y-^M C0S9)| + 


I + i 


(y-i8M^cos^0)^ 
izVa^: g2i(y-/9M^os^0) 


I _g2i(y-iSM^cos^0) 


2E 




(y-)8 ^coidf 


(y-^M^OS^) M^CO^d 


^ 2 ( y+)8 M^cos^ 0) ][ 


(39) 


(to) 


,4iy 


2v^ e 


2i(y~/3M^cos^^) 


+ 2 ^ — — — 

"rr y-^^M^cos^d ( y-yS M^cos^0 ) \/y\(B M^cos% 


2(^+jSm2co^0) 

+ Tl 


Equations ( 36 ), -(39), and (4o) can be used to determine the aerodynamic 
response function for large values of 7^. 


Limiting Behavior as |y^| -»-oo 

The limiting form attained by Eq. (39) for very Jarge values of l7^| 
can be derived by expanding Eq. (39) in powers of l/y2. The result is 


Tl = 


j_ 

7T 


_J 

^{/B+y) {y-/3M^coie) 


(to) 
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provided 0J^cos^d ^ 0. Equation (4l) shows that is proportional to 

since hoth 7 and )8 are linear functions of K. However, if 7 - fil^cos^d¥ 0 
Eq* (39) will approach 


T. 2(l + i) 

*■ 7ry7r(^+x) 


(h2) 


1 

for large values of 7. This resTilt is proportional to K ^and occurs when 
the gust phase velocity along the span is sonic relative to the mean span flow 
velocity Usin0. This phenomenon will occur when the wave front encounter 
angle is equal to 


= cos“' 




1 + 


M^OS% 


0< M < I 


(43) 


A comparable expansion of Eq. (40) for very large values of 1 7^ I yields 
an expression for which is proportional to K”3/2^ provided 7 - /3^^cos^0 * 0. 
However if 7- pi^coe^d = 0, t]^e limiting form attained by the expansion of 
Eq. (4o) is proportional to K" 2 . This large variation in the asymptotic 
limit of the moment transfer function can, for example, have an effect on the 
unsteady aerodynamic response of a helicopter rotor blade encountering the tip 
vortex shed by a previous blade. 
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DISCUSSION OP THE APPROXIMATE SOLUTIONS 




The ohjective of this section is to demonstrate the validity of the 
approximate solutions that were developed in the previous sections for the 
aerodynamic response functions. First, consider the eqiiations for the lift 
transfer function, Tj,, for both small and large positive values of y (i.e., 
Eqs. (27) and (39) t respectively) and compare the two-dimensional limits of 
these equations to the exact results of Refs. 5 and 6. This comparison is 
presented in Figs . 4 and 5 in which the complex function Tl is Pioi-ted in 
the phase plane- at specific values of nondimens ional wave mamher, hK, for . 
Mach numbers of O.3 and 0.6, respectively. The solid curve appearing. in 
both figures corresponds to the exact results presented in Ref s . 5 and 6, 
while the dashed curves represent values computed from Eqs. (27) and. (39) 
and. a linear curve fit procedure to Join the upper limit.. of validity of 
Eq. (27) (i.e., 7 - 0.2) to the lower limit of validity of Eq. (39) .(i-e., 
y = 0.7). Figures U and _5 show that the present approximate solutions for 
the lift transfer functions are in good agreement with the exact result 
over a wide range of bK. The accuracy of > Eqs. (28) and (40) in predicting 
the two-dimensional moment transfer function is shown in Figs . 6 and 7 for 
Mach numbers of O.3 and 0.6. The values of the moment transfer function 
for the region between the upper limit of Eq. (28) (i.e., 7 = O.l) and the 
lower limit of Eq. (ko) (i.e., 7 = O.7) were determined by employing the 
same linear curve fitting procedure that was used for the lift transfer 
function. The agreement between the approximate solutions and the exact 
solutions for the moment transfer functions is seen to be comparable to that 
attained for the lift transfer fvinctions . 

The accuracy of the approximate expression for the lift and moment 
transfer functions governing the region and zero angle of sweep is 

first demonstrated at = 0 by comiiaring the results predicted by Eqs . 

(34) and (35) at this limit with those predicted by Eqs. (27) and (28). In 
Eqs. (27) and (28) the limit as y-*0'*' yields exact results while the limit 
of Eqs. (34) and (35) is approximate. A comparison for y^ = 0 is shown in 
Figs . 8 through 11 in which the complex functions T^ and Tm are plotted in 
the phase plane for specific values of nondimens ional gust wave number, bK. 
Once again the agreement between the exact solutions and the approximate 
solutions is good. To establish the accviracy of Eqs. (34) and (35) for 
values of 7^< 0 a comparison was made with the results obtained from an 
analysis based on the theory of Ref. 7. This comparison is presented for 
the real part of T^ and Tm in Figs . 12 through I5 for an encounter angle 
of 90 deg, zero angle of sweep and for Mach numbers of O.3 and 0.6. (The 
imaginary part of T^ and Tjjj predicted by Eqs. (34) and (35) is zero, which 
is the correct solution for this particular encounter.) The solid curves 
appearing in these figures correspond to the exact resvilts while the dashed 
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curves were predicted by Eqs . (3^) and (35 )• In addition, the resiilts 
predicted hy Eqs. (39) and (4o) are included as the dot-dash curves. The 
agreement between the approximate solutions and the exact solution is very 
good for values of i^>1.0. For values of Kb s 1.0 the agreement between 
Eqs. ( 34 ) and (35) and the exact results has deteriorated scsnewhat but is 
more than adeq\iate for computational pijrposes . 

’The validity of these approximate solutions is clearly demonstrated by 
the good agreement with exact theory, and the closed form structure of the 
solutions facilitates their use. The latter is particiilarly important when 
these results are applied to specific problems of unsteady aerodynamic re- 
sponse where computer time and core storage may be limiting factors . 
Furthermore, the mathematical procedure that was employed to develop the 
approximate expressions for the unsteady response functions is sufficiently 
general to be applied to additional unsteady aerodynamic problems . In 
particular the current procedure can be used to predict the unsteady- forces 
generated by the pitching or plunging motion of a two-dimensional airfoil 
in a subsonic stream. The derivation and discussion of these transfer 
functions is presented in Appendix IV. 



APPLICATION OF THE APPROXIiyiATE .SOLOTICMS , . 

The theory of the preceding sections has been, .applied. to .the study of 
the lift and moment response of a helicopter rotor blade .encountering the 
trailing vortex from a previous blade. The model that has. been developed, 
to simulate this aerodynamic interaction problem assumes the rotor blade to 
be an infinite swept airfoil. The trailing vortex is modeled as a recti- 
: linear vortex. translating at a velocity e,qual to the. mean stream velocity. 

The lift and moment coefficients per unit span generated by the encounter 
are obtained by the Fourier superposition of the lift and moment responses to 
an infinite series of oblique gusts which. constitute the. Fourier decomposition 
of the rectilinear vortex.. The mathematical, details of this procedure ""are 
discvissed in Appendix V, while the computer program used^. to compute the " 
example cases is described in Appendix VI. A series of computations were 
performed to determine the effects of , encounter angle, Mach; number, and 
vertical vortex displacement on the lift and moment, response , and- these 
results are presented- in Figs. l6 through .23. In all cases examined, ,.;bhe. 
airfoil sweep angle was fixed at zero (i.e*, the incoming flow was.,normal, 
to the leading edge line of the wing). ,. The lift and moment coefficients 
per unit span appearing in these figures- are normalized, with respect to the 
parameter r/Ub where F is the circulation of the free vortex,.. IJ is. the free 
stream, velocity, and b is. the semichord of the airfoil. , A tyjpical range of. 
this parameter for a helicopter blade is. 0.2 sr/Ub <2. -.Figure- 16, shovs. the. 
lift coefficient per unit span generated- by a parallel encounter, of a. yortex 
with a two-dimensional wing as a function of free stream tfech nmber and , 
nondimens ional time, t’ = Ut/b. At t* = 0 the position of the vortex is 
one semichord below the airfoil leading edge, while at -t-' ,= 2. it is located 
one semichord below the trailing edge..- Ibe lift coefficient calc.uia^:,ed on 
the basis of . quasi-steady lifting line .theory, is also shown as , a, dashed, line 
in Pig. 16. This- result is presented; ^^;o illustrate the effects of. .unsteadi- 
ness . It is seen from the figtire that lifting line theory overestimates , the 
maximum values of the. lift coefficient when. the vortex Is directly beneath 
(0^t':^2) or slightly downstream (t'> 2) of the airfoil. This behavior can 
be attributed to rapid time variations, in downwash that -occur during the time 
the vortex is beneath .the airfoil . In Fig. IT the moment coefficient about 
the quarter-chord produced by a parallel encounter is shoyn .as a function of 
Mb ch number and nondimens ional time t'.., A non- zero moment coefficient about 
the qviarter- chord for .a parallel encounter can be directly attributed to 
compressibility, since the moment is zero for an, incompressible' fluid.' (This 
is shown on page 28 t of Ref. 9« In addition, it should be, noted that lifting 
line theory - does not predict, pitching .moment . ) It , is seen, from Fig.' 17 that 
significant values of (r/Ub) can be obtained, which increase with 

increasing Mach number. 
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The results appearing in Figs . 18 and 19 illvistrate the effects of 
three-dimensionality at- a free stream approach Mach nmher of 0.6. These 
figures represent the lift and moment coefficients, respectively, produced 
hy a vortex encountering ah -airfoil at an oblique angle'. as shown in the 
sketch ih Pig. l8. Both resixlts are presented as a ftmction of nondimen- 
sional time t', 'which is more generally defined as 

' “ ' t' = ^ cos (0-a) sin a 

D D 


Ih this equation a is the angle between the vortex and the leading edge of 
the wihg,and Xo is the spanwise coordinate. The spanwise position of the 
vortex at time t = 0 is x^ = 0, while its vertical position is one semichord- 
below theVing. The lift response shown in Fig. l8 appears to be a weak 
functioh'"of the encounter angle, while the corresponding moment response, 
shown in. Fig. 19, appears to be very sensiti-ve to variations in the encounter- 
angle. This effect can be directly related to the asymptotic behavior of , 
the Fourier transform of the moment response . Over most ■ of the range of 
encounter angle, a , from parallel to normal ehcounter, this asymptotic limit 
is inversely prbportiohal to ■ (details can be found in the section 
entitled Limi-ting Behavior as ly^j -*•<») . However, if the value of the 
encounter angle is^such that a = a* (cf. Eq. (43)) the asymptotic limit is 
proportional to K"2 . This was the case for a = 30 deg. It appears from 
Fig. 19 "that there is a definite transition in the' character of the moment 
response' at a = a*. ' 

A number of computations were performed to determine the effect of the 
vertical position of the free vortex beneath the airfoil on the unsteady lift- ' 
and momeh't response generated by an encounter. The resiilts for the normalized 
lift and’ moment coefficients per unit span generated by a parallel encounter 
at a Mach mimlier 0.6 are shown in Figs. ,20 and 21 as functions of nondimen- 
sional tiiK t' and nondimens ional vertical distance h' = h/b (where h is the 
vertical height of the free vortex beneath the wing). The geometry of the 
encounter at time t = 0 is shown in the sketch of Fig . 20 . These results 
indicate tha't the effect of increasing the vertical distance between the 
airfoil and' the free vortex is to -reduce the lift and moment response . Fig-^ 
ures 22 and 23 show the maximum positive and negat I've values of the response 
curves as functions of- vertical displacement and encounter angle. These • 
results were obtained for an approach Mach number 0.6. Decreasing the 
vertical displacement between -the free vortex and the wing causes the lift 
and moment response to increase for all encounter angles examined to date. 
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These results also show that the maximum positive values of the lift and 
moment response curves appear to he relatively insensitive to the encounter 
angle. However, the maximum negative values are strongly dependent on the 
encounter angle . 

The final set of results that are presented is for a UH/1D rotor in 
forward flight at 90 Kts where the advancing blade in the second quadrant- 
encounters the tip vortex shed by the preceding blade. Data for an ex- 
perimentally meastired wake geometry (Ref. 12) has been used to predict, the 
spanwise distribution of the aerodynamic loads when the vortex-rotor blade 
intersection point is at O. 85 , 0.753^ and O.Sjk of the rotor radius. The 
vertical displacement of the vortex beneath the plane of the rotor is 0 . 87 ^ 
semichords while the azimuth position of the Intersection at O .85 of the 
rotor radi\is is 92 deg, at 0.753 of the rotor radius is 104 deg, and 
at 0 . 67 ^ of the rotor radius is II 5 deg. The resulting normalized lift 
and pitching moment coefficients are presented in Figs . 2h and 25 for the 
three blade radial locations. The results indicate the lift and moment 
coefficients to be relatively insensitive to radial location. This result 
can be attributed to the small variation of the geometry of the encounter 
with radial position. A variety of blade/vortex Intersections can occur in 
a rotor depending upon advance ratio, azimuth position, lift and propulsive 
forces, and maneuvering condition; however, an eval\iatlon of such intersec- 
tions is beyond the scope of this study. 


26 



RESULTS AND CONCLUSIONS' 


A primary objective in developing the present analysis was to sig- 
nificantly reduce the computational time required "by the theory of. Ref. 7 
to predict the aerodynamic response resulting from the encounter of a two- 
dimensional swept airfoil' with a three-dimensional oblique sinusoidal gust 
in; a compressible 'stream. This general objective was achieved; and the ' 
folldwing are specific results obtained and conclusions, reached in the 
development of this analysis; 

, r. A simple set of approximate expressions was derived to model the- 
airfoil unsteady pressure, lift, and moment response. ■ 

2 . In predicting the airfoil; aerodynamic response, use. of , these 

- approximate expressions significantly reduced the computation time 

• (by a factor of 10 or more) over that required-by the theory- of Ref. 7* 

. -'3. The validity of these expressions was demonstrated by the close 

- agreement with; the exact results, (i .e based on the , theories of Refs. 5 
- ;-and 7);- over a . wide- range of values of the. governing parameters ,. 

U. The approximate expressions clearly show the influence of compres- 
sibility and three-dimensionality on the response functions, and they 
show that a major effect is to generally increase the asymptotic decay 
rate of the response functions relative to their incompressible two- 
dimensional counterparts . 

The approximate expressions were used to predict the aerodynamic 
loading on a helicopter rotor blade encountering the tip vortex shed by a 
previous blade. The major conclusion reached in this application of the 
analysis is that compressibility, unsteadiness, and three-dimensionality are 
important terms which must be included to correctly predict the aerodynamic 
response, particularly if the vortex is located within two chord lengths of 
the airfoil. Some of the specific results that were obtained in this 
application are: 

1. The influence of unsteadiness is to reduce the magnitude of the 
lift response compared with that predicted by lifting line theory. 

2. The influence of compressibility and three-dimensionality is to 
cavise the generation of a moment about the qviarter- chord which would 
be zero if the flow field were incompressible and two-dimensional. 
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3 • The lift and moment response amplitudes increase as the Mach numher 
of the flow increases . 

4. The maximum value of the lift and moment response decreases,, rapidly 
as the distance between the center of the vortex and the plane of the 
airfoil increases . ' ' , 1 



RECCMIENDATiaNS 


The tasks described in the following recommendations will provide a 
framework for futiire studies of the vortex encoianter problem: 

1. An exijerimental program should be undertaken to measure the 
aerodynamic forces and moments generated during the encounter of a rotor 
blade and a vortex. These measurements should be made over a wide 
range of flow conditions and specifically should include the region of 
vortex- induced stalled flow on the suction surface of the rotor blade. 

2. The present analytical model does not account for the influence 

of a tip region on the forces and moments generated during an encounter. 
Work, should be undertaken to introduce this effect into the present 
analysis . 

3» The results predicted by the present analysis Indicate that lifting- 
line theory is,, inadeqiiate for predicting the aerodynamic forces and 
moments generated during an encounter . It is therefore suggested that 
rotor blade dynamic response analyses which rely on lifting- line models 
to determine the aerodynamic forces be modified to incorporate the 
present analysis to yield a more realistic prediction of the forces and 
moments generated during an encounter. 

4. Finally, the present analysis should be combined with an aco\istic 
analysis to predict the acoustic field generated by a vortex encounter. 
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APPENDIX I 


APPROXIMATE EXPRESSIONS FOR THE AERODYNAMIC 
RESPONSE FUNCTIONS FOR SMALL POSITIVE VALUES OF 7 


In this Appendix approximate expressions for hoth the circulatory and 
noncirculatory solutions of Eq. (l6) are derived for small positive values 
of 7. Consider first the variable associated with the noncirculatory 

flow field which is obtained by solving Eq. (l6) subject to the boundary 
conditions, Eqs . (l9) and (20 ). From the results presented in Ref. 7 it 
can be shown that a solution for the noncirculatory flow field of order 7 
(i.e., a series which neglects all terms of order 7^^ for n>l) is given by 
an expression which satisfies Laplace's equation. 

v*4.‘^’=0 (ll5) 


and the boundary condition, Eq. (19)* The solution of Eq. (U^) is given in 
Ref. 7 as 

=y^ Ansinrvue"^^ (l<.6) 

n=i 


where the elliptical cylindrical coordinates tj and f are defined by the 
equations ' 


X, = cosh Ceos 17 


C^7) 


and 


X2= sinh C sini7 


(48) 
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The surfaces ^ = constant and ■»? = constant are elliptic and hyperbolic 
cylinders, respectively, and the degenerate ellipse i = 0 corresponds to 
the portion of the plane occupied by the airfoilj that is, -IsXq^sI. 

The constants in Eq. (46) are determined by requiring to satisfy 

the botindary conditions in Eq. (19)- 


dX 


(I) 


2 


d4> 


(I) 


X2=0 




l_ 

sim7 


-E 

n - I 


ApO sin T)V) 
sin‘>7 


Eqs • ( 19 ) and 


(49) may be rearranged to yield 


00 

^ Apnsin T\r) 

nz\ 


b 



M^COS^0 


sin'J^Uge 


iiScos^? 


(49) 


(50) 


If both sides of this eq\jation are multiplied by sin m’>7 and integrated over 
the interval 0 :£ '>?:S7r the result is 


bu. 


I 

7T 


A - 

- 7t 


bu. 


^ I - 0 


. i)8cos^ 

J sinT7 Sin m‘»7 e 


/ ■’'■f I ificos^7 

I cos(m + i)'»7 + cos(m-i)‘J7j e d**? 


(51) 


^ I - M^COS^S 0 


Upon introducing the integral relationship 


m > I 


Jp{;9) 



locos'*? 


cos nT7 d*i7 


(52) 


and the Bessel Function recurrence fonnula 


Jp.,(Z) + Jp.,(Z) 


2n 

z 


Jn(Z) 


(53) 


31 



into Eq. (5l)> one obtains 


' M 


2bUg 

I - M^COS^0 


. n- I 


/? 


Jm (/3) 

M 


m > I 


(5'+) 


which when combined with Eq. (^6), transforms Eq. (46) to 


<J> 


(I) 


2 bU2 




4- Sin 0 77 6'"^ 


(55) 


. An expression for the noncirculatory pressure field on the surface <of 

the airfoil is obtained by combining Eqs . (2l) and (55) yielpTing i ■ 

/ 


S(l). I 

* ’ (i-ivi^cos^^ sin 17 


cos 17 +i J,()3)j X 


i cos 17 + I cos (6-o)tj 


( 56 ) 


Since sin77 = 0 at the leading and trailing edges of the airfoil, the pressure 
distribution given by Eq. ( 56 ) is singular at these points. 

• r '• ■ ■ ; . , ' 

The trailing edge singularity is removed by adding a circulatory field. 
The approximate expression for the circulatory flow field pressirre distribu- 
tion on the surface of the airfoil for small positive values of 7 is obtained 
by expanding the exact solution for the- circulatory pressure distribution 
presented in Ref. 7 in a power series of 7. This pressure distribution is 
equal to 


8(2) -(2) i[-^M^COS^0COST7^-jrj54^ >^3-KUCos(^^^^ (57) 

P = P e 
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is defined "by Eq. (62) of Ref. 7 as 


where 



(0,17) = - 


PqUcos d 


conj 


CO 


£ Bncenln.q) (0,q) 

n=0 


q=y% (58) 


and the symbol conj denotes the complex conj-ugate operator. The coefficients 
which appear in Eq. ( 58 ) are defined hy Eq. (5^) of Ref. 7- For small 
positive values of 7, an expression for p^^^> correct to first order in 7, 
can he obtained by retaining the first two terms of the series in Eq. ( 58 ), 
namely 


X 

(59) 


■ ^ , yre, Me,‘2>'(0,q) 

Jz ceQ{r],q) + — rgp- — -ce(77,q) 

®o Meo ' (0,q) 


When ceo iVfO.) and cei (’?>q) are expanded in a power series, the following 
resialts are obtained to first order in 7 : 

~ ( 60 ) 





/OqUCOS© 

b sin T) 


conj 


■/T 


(0,q) 


ce, (77,q)=cos77 (61) 

yFB| 

Similarly, if the quantity Bq expanded (i.e., by employir(g the 

results presented on pages l46 and 147 of Ref . 7 with 7 and 0 replacing ju, and 
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Vf respectively) the result will be 


( 62 ) 


/{ 2 ) 
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The first order expansions of the remaining Mathieu f met ion are; 
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( 66 ) 
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Upon substituting Eq. (63) through (66) Into Eq. (62) and neglecting terms of 
higher order than 7 which appear in the denominator of Eq. (62) the result is 




(y/2)-i(rM^)Hj^^(i9) 


(67) . 


Equation (6?) is reduced to its final form by interpreting 1 + iyS as the 
.leading terms in the expansion. of e^^. With this replacement; Eq. (67) 
becomes 


yTB,Me,'<2) ^ i + i r^/0 ^ / 2 ) - 

B(^eo'<2) ” H|(2)(|S) 


(68) 


The substitution of for (l + iyS) causes Eq. (68) to approach zero as 0-*oo. 
This substitution improves the accuracy of the approximate expressions for 
the lift and moment transfer, function for large values, of >3. 

Equation (68) contains an expression of apparent second order in 7. 
However; if this term is neglected; a term of first order in K would be lost. 
This first order term results from a combination of the expressions within 
the brackets and the Hankel function (/3). Thus if 7^ is neglected in 

Eq. (62) the validity of the approximation is reduced to zero order in K. 

When Eqs . (60); (61); and (68) are introduced into Eq. (59) an approximate 
expression for p'^' is obtained 


bsinr) 


I-H- 


h/2)(/3) 


-COS17I 


(69) 


35 



where the constant C must he assigned a value which will cause at the 

trailing edge to he equal in magnitude hut opposite in sign to evaluated 

at the trailing edge . 

The final form for the pressure distribution on the surface of the airfoil 
for small positive values of y is obtained by substituting Eq. (69) into Eq. 
(57) and adding the result to Eq. (56) to yield 




^pUtTgCOS e 

V l+X| 


Jo(/ 3 )-iZJ,(/ 3 )‘ 


+ z 


(TO) 


-(Sn^co^d X| 4 K sin ax^- [l? u cos {0-a) t|j 


where the variable 


Z= - 


i [hq” (y/ 2 ) - Ho^'H/ 3 )] 


( 71 ) 


The equations for the lift and moment transfer functions T^ and Tj^^, 
respectively, can be obtained by integrating the zero and first order moments 
of the pressure distribution over the chord of the airfoil. The resulting 
expression for the lift transfer function T^^ is ; 


T, = e 


ibKcosa 


Jo(;g)-iZJ,( 6 g)] 


If Z 


[jq( m ^/3 cos^e) f- i J , ( M^/3 CO ) ] (72 ) 
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while that for the moment transfer function Tjyi about the quarter-chord is: 


- ibKcosa r 1 

^^>1/4'-® fTz j"" 

[~Jo(M^/3cos^0)+ J2(WI^/3 cos^ 0) - 2i J|(M^j8cos^e)] + Tl 


(73) 


The factor which appears in Eqs . (T2) and (T3) bas been added to 
shift the phase of the lift and moment response relative to the gust field 
to the leading edge of the airfoil. 
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APPENDIX II 


APPROXIMATE EXPRESSIONS FOR THE AERODYNAMIC 

..... 2 

RESPONSE FUNCTIONS FOR NEGATIVE VALUES OF 7 


The ohjective of this Appendix is to present the derivation of the ' 
approximate expressions for the aerodynamic response functions for the 
region T^sO. This derivation was accomplished "by employing an extended 
version of the similarity rules developed hy Graham (Ref. 5); which were 
modified to include the effect of airfoil sweep angle. Graham shows that - 
.the solution' for the modified velocity potential 4> (cf . Eq. (15))' for the" 
region 7 0 and for zero sweep angle is directly related to the solution ' ■ 
for ♦ associated with the prohleiji of an unswept two-dimensional airfoil en- 
countering a three-dimensional oblique sinusoidal gust in an incompressible 
stream. This correspondence also exists for swept. infinite airfoils and is 
expressed by the equation,;. 




y / 1 - 


(74) 


which is derivable from Eq. (l6). Here ♦('y,/3)njC solution to the 

three-dimensional incompressible problem expressed in terms of y and/8. 

The substitution of Eqs . (t 4) and (l3) into Eq. (2l) yields the. expression • 
for the compressible pressure field 

/ c . 

P /“ 5 X . 

b I- M COS 0 

exp [-i/3^M^cos^0X I +iKsinaxj- i Kucos(0-a)t] 

The equation for the corresponding incompressible pressure field is obtained 
by letting M-*0 in Eq. (75) , ■ 


d4> 


INC 


dXi 




(75) 


pir,J3) 


p^ucose 


INC 


exp 


■aTT 


(76) 


iKsin 


0X3- iKUcos 


(e-a)t] 
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and it is seen that Eqs . (75) and (t6) are related "by the equation 


P^ = 


P(y.)S)iNc 

V^i-M^cos^0 


(IT) 


which is a pressure similarity rule comparable to Eq. (7^) for the modified 
potential. 

A formula for the pressure field p(*>'^y9)™c derived from the 

approximate expression developed by Filotas (Ref. 4) for the pressure 
distribution on the siirface of an unswept airfoil encountering a three- 
dimensional gust in an incompressible stream, 


(I) 


(2) 

g-KpC0S/3p 


Iq( KpCOS /Sp) 4*T|( KpCOS /8p) 

( 3 ) ( 4 ) 

{exp (iKpC0S/3pX3)}x (exp{- iUKpsin/Spt} 


I -X, 


i + x 


(78) 


where 


-I <X,< I 


H(Kp,3p) “ 


j 

[ / ^ 



” iKcj 

e '^1 

f 

|sin/3p -7T)8py + cos — 

\ / /^ \ 
j i+27rKpy-fcos-g^| 



1 + 7T Kp( 1 + sin^/Sp t 7T KpC0S)3p) 


(79) 


The subscript F which appears in these equations refers to the variables 
defined in Ref. 4. The transformation of Eq. (j8) to Pjjjc is accomplished 
by letting the terms in brackets (l), (3), and (4) become 

bracket (l) 




( 80 ) 
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"bracket (3) 


{exp[iKFCos^rxJ| — {expfiKsina xj} 


( 81 ) 


"bracket (h) 

{exp[-iKpSin/3pt]} — {exp[-iRucos(0-a)t]} ( 82 ) 

while the remaining terms which appear in "bracket (2) are transformed "by- 
replacing the varia"bles Kp and ySj. "by the expressions 



(83) 


= ir/Z -8 (84) 

where 

8=ton"'y/3 (85) 


Introducing these transformations into Eq. (78) and combining the result with 
EQ* (TT) yields an expression for p® 


yD^Un^COS^ [ 

exp [-iKsinaX 5 -iKUCos(e-o)t] e 


»’yxi 


/I -X 
V i + x 


i- X 


2 2 
COS ®X| 


( 86 ) 


where 


H(y./3) = 


. /-T i 6 ir[7r/2-dfl+8in8/d \ 

e r* * Li*.ln8/2] f 

,^i + Tyy 2 +jQ 2 [ii-cos S+TTV-ye+^isins] 


(87) 


\ 


ko 



The lift and moment transfer functions per unit span are obtained hy 
integrating the appropriate moments of the pressure distribution, Eq. (86), 
.oyer the chord of the airfoil. The expression for the lift transfer function 
Tl is 


Tl 


Io(/ + i^lVI^COS^0) + COS 0) 

ioir) + I.ty) . . 




( 88 ) 


while the corresponding expression for the moment transfer function about 
the quarter-chord is 


^ „ iKbcosa 

S/4=~2e 

+ ii (y+ i /3 m^cos^0 ) 


h(x,/3) 


IJ[y+i;3M^cos^0) 


y + iiSM^cos^0 ) 


+ Tl 


(89) 


The factor has been introduced to refer the phase of the lift and 

moment response fxmctions to the gust downwash at the leading edge . 



APPENDIX III 


SOLUTION FOR LARGE VALUES OF |7^ 


In the main "body of the text it is suggested that for large positive 
values of 7^ the pressure distribution on the surface of the airfoil is 
composed of two cylindrical acoustical waves whose origins are the leading 
and trailing edges of the airfoil. Furthermore, it is stiggested that the 
solution for the pressure field whose origin is the leading edge could be 
approximated by the solution for the pressure field resulting from the 
encounter of a three-dimensional gust with a semi^-inf inite plate. The 
solution for the modified velocity potential 4* associated with this pressure 
field satisfies Eq. (l6) and the boundary condition of Eq. (l9) on the surface 
of the plate. , The solution for <i> can be obtained from an analysis based on 
the Weiner-Hopf technique, the details of which can be found on pages U8 to 
98 of Ref. 13- This solution is 


*(X,,0) 


E(iti)e'^" 


b U 2 

■/ 1 -M^COS^ - 


X,'>0 


( 90 ) 


where 


^ ^ ^ [^{p-y)'K \ ]“ ' ^ ] 


( 91 ) 


and C and S are cosine and sine Fresnel Integrals, respectively. The pressure 
distribution is obtained by substituting Eq. (90 ) into Eq. (21 ) yielding the 
result 


^UUgCos^ (i + i)exp 

p(l) = — . " " ~ T~ ^ 

V I - M^cos^S 

[i(X-M^/3cos^0)x'+ iKsinaxs -iKUCos(0-(j)t] 


(92) 



where x{ = X]_ + 1 and the superscript (l) denotes that Eq. (92) Is an 
approximate expression for p® in a region removed from the trailing edge 
of the airfoil. The magnitude of the pressure dlstrlhutlon as given hy 
Eq. (92) monotonically approaches zero as the distance from the leading 
edge of the airfoil "becomes large without "bound . However, the value for 
the pressure dlstrl"butlon predlctedjDy the exact theory of Ref. 7 for the 
region downstream of the airfoil (l.e., X^^l, X2 = O) is zero. Hence, 
the pressure dlstri"bution, Eq. (92), is incorrect in a region surrounding 
the trailing e(^e of the airfoil, the magnitude of the error "being propor- 
tional to This error may be corrected by adding to p'^' a pressure 

field p(2) which is equal and opposite to p'^/ downstream of ’ the airfoil 
(i.e., X^^al, X2 = 0) and has its normal gradient vanishing everywhere on the 
airfoil. An approximate expression for this press\ire field can be developed 
for'y^> 0 "by considering the problem of an oblique cylindrical acoustical 
wave impinging on an infinitely compliant, semi-infinite surface (i.e., a- • 
surface on which the pressure vanishes) in a uniform compressible stream. 

(Por y^< 0 the problem of an oblique cylindrical cut-off wave impinging on 
a compliant surface can be , considered . ) The solution for the scattered 
pressure field p'^' is developed by again employing the Welner-Hopf technique 
(Refi'13). The construction of the solution for the pressure field p'^^ is 
accomplished by Introducing a modified acceleration potential H which is 
defined In termsof p'^/ by the equation 



i [-/jM^cos^© Xj + Rsinax3“KUcos(0-a)t] 
yoHe • 


( 93 ) 


where the variable transformation X^ = X-[ -.2 now places the local origin at 
the trailing edge. Substitution yields the result that H satisfies Eq. (I6) 
everywhere in the flow field. 

Downstream of the airfoil (l.e., X^aO, Xg = O) p^^^ is equal in 
magnitude but opposite in sign to p(i). Hence from Eqs . (92) and ( 93 ) H must 
equal 


H = - 


UugCOSfi 


(l+i) e 


iy(x(+ 2 ) 


-2i)8M^CpS^0 


^ I - M^cos^fi ^ 27t(X| + 2)(^+ r) 


(9U) 


>0, X 2 = 0 


k3 



(95) 


£ownstream of the' trailing edge. Upstream of the trailing edge (i.e., 
0, X 2 = O) the normal gradient of p(2) must vanish, hence 


an 

ax, 


= 0 


X, < 0, X, = 0 


The solution for the pressure field p^^^ may he derived from the solution 
of the mixed boundary value problem for Hj_ which can be found on pg. 79 of 
Ref. 13 . The value of H on the boundary 0, X 2 = 0 is given by the. 
expression 


H 


3 iir { CO 
I - 3/2 - 

'"O 


r® _/ -iru r d r°°-i/2 ix(u*f) 

Jm(u, x„ 0 )e du|-^jTe e f(U + C)dC 


3ii 

I -3/2 ~z~ I r'*' -I7U d 
-_±rr / Me -|j-I(U)dU 


00 _ 


iru cJ 


(96) 


where 


X| < 0, Xp = 0 




M(U, X,, 0) = 


(U -X.) > 0 


(U-X,) < 0 


(97) 


and the variable f is equal to H evaluated on the boundary X, > 0, Xp = 0 



The evaliiation of the inner integral in Eq. (96) is accomplished hy first 
introducing Eq. ( 98 ) into Eq. (96) and replacing H(U + withEq. ( 9 ^)« 
These manipulations yield 


I(U) = - - 


Uu,cos9 


2i(y-)(3M^cos e ) 2iy(u+£) 

(I + i)e j"® e ^ 


a/ 2 tt ( / 3 +y ) 0 y/u + C+ 2 ‘sy~f 


- dC ( 99 ) 


This result may he reduced to 


I(U) = - 


UOgCOs 9 


. 2i( y- ySwi^cos^0) 
le 


2iyu 


^\ - V2y{(3+r) Vu + 2 


(100) 


hy employing the method of stationary phase (i.e.,'Ref. li<-) in evaluating the 
integral over Substituting Eq. (lOO) into Eq. ( 96 ) and performing the , 
differentiation with respect to U yields 


H = 


UOg COS 9 
^ 1 - M^COS^S 


. 2\{Y- Sm^cqs^B) 
le 

VzrifB + r) 


(101) 


. ^ 2iyu . 

g-iyxifoo e 

2iy- * 

*^0 V(U-x|){U + 2) 

2(U + 2) 


dU 


x( < 0 



2UQ,cos9 


H = 


Vi-«'cos^e 


( 102 ) 


iyx, /.CO 


2 iyT‘ 




172 


2i7- 


2(t^ + x'+2) 


dr 


1+5 



An approximate value of the Integral is obtained by noting that for large 
values of y the major contribution of the Integrand to the integral comes 
from the region surrounding the lower limit* Thus, Eq. (102) may be 
approximated as 


UU9C0Sd 

H “ . - , .. - 

"^/l -M^COS^0 


, 2i(y-M^/9COS^0) . 

le piyx'i 

yy(/3 + y) 



UQgCOS© gjy(2+^|) ■ -2iM^ygcoS^g 

,yi-M^cos^0 y7T(/3 + y) 



The pressure distribution p(^) on the boundary 0, Xg .= 0 is obtained by 
combining Eqs. (93) and ( 103 .) which yields 


( 2 ) 

P 


-^UugCOsS 

*v/l-M^COS^0 


i{[y-/3M2 cos^s] xJ+KbsinXj-UK cos(5-a)tj' 

^ ~r==^=- ^ 

.yr(^Tyr_:,\ 


(io4) 


-^-e(v'2)'(2-x;))J osx;<2 

An approximate expression for the pressure field p® is obtained by combining 
the expressions for p'^', Eq. (92) and p'^', Eq. (104). The result is 


P* = 


/ouagcos© 


I + i I + i 


yi-iviWa ■v/7r(/3+y) (y 2 X,' 


+ E 




X ( 105 ) 


y-/9M^cos^9 


x| + Rb sin a X,- UKcos(9-a)t, 


U6 


exp I 



In Ref. 7 it is shewn that upstream of the airfoil the solution for 
p® must vanish. However, the approximate expression given hy Eq. (IO 5 ) 
equals p® = p^^^ upstream of the airfoil since the pressure distribution 
p^^' vanishes on the boundary 0, X 2 = 0. It is seen from Eq. (10^4-) 
for p(2) that the magnitude of this expression for large values of 
is proportional to 7"^. Hence^ although the pressure distribution, 

Eq. (105) is correct near the trailing edge, the result p® = shows 

that it is now incorrect in the region upstream of the airfoil leading 
(i.e., X]^< -l), the error being proportional to 7 "^. This error 
could be removed by employing an analytical procedure which is identical 
to the one used earlier to remove the error associated with the 
pressure distribution p^^^. The result wo\ild be an estimate for the pressure 
distribution p® which would now be in error in the region surround iM the 
trailing edge, the magnitude of the error being proportional to 7“^' • This 

could be further corrected at the expense of the leading edge solution, and 
if desired the analytical procedure could be employed as many times as 
necessary to achieve the desired degree of accuracy. The final form for the 
pressure distribution p® on the surface of the airfoil wo\old then be an 
asymptotic series involving powers of 7“ 1 / 2 . In the present analysis the 
series is terminated after the second term, p'^'', which yields, the result 
given by Eq. (105). 

The approximate expressions for the lift and moment transfer functions 
are obtained by integrating the appropriate moments of the pressure distribu- 
tion, Eq. ( 105 )* over the chord of the airfoil. The resulting integration 
yields an expression for the lift transfer function Tj^ which is 


Tl = 


y 7 T() 3 +y) 


(i + i) 


M^cos^ 


E [«/ 2 ] 


+ [, _g2i(y-^M^C0S^e) J 


I- 1 


2(/-^m2cos^) 


(106) 
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[ysy] 


i yir 


{y- 


2i(y-/3i^col0) 

_e 

y^r^SM^cos^ 


[y 2(X+-/3 m2coS?S) ] 


UT 



while the corresponding expression for the moment transfer function about the 
quarter-chord is 






'^^(y-y3M^cos^9) 


I +i 


{y-pM^cos^e)^ 


I - e 


2i(y-/9M^COS^0) 


" 2E[2yy"] 


\zyTF e' 


2i(y-/3M^c6s^0) 


(y-/9MW0)^ (y-/9iviW0) yrt/sM^cos^e 


X 


[V2(y+/3M^COS^5) 


2 - 


2{r+>9M^cos^0) 


(107) 


+ 2 


2^11 

i V tt 


4iy 


2yrr e 


2i(y-)8M^COS^0) 


r^/3 M cos 0 (r-^M^cos^0)^-v^y+ /3 m^cos^0 


>y^2(y+j9M^cos^0) 


+.T, 


The corresponding expressions for the lift and moment transfer fxmctions for 
y2 < 0 are obtained by replacing y by iy in Eqs . . (106) and (lO?)* 



APPEnroix IV. 


TWO-DIMENSIONAL PITCHING AND PLUNGING MOTION OF 
/■. AN AI^OiL .m' A COMPRESSIBLE STREAM 

The mathematical procedure developed in the previous sections was shown 
to he capable of generating simple approximate expressions for the aerodynamic 
response functions. .associated with the gust encounter problem. This result 
might lead one*, to attempt to hpply the same procedure to develop aerodynamic 
response functions for arbitrary motion of a finite lifting siorface . This 
general three-d imensional; problem , cannot be analyzed exactly by the present 
theory because tip effects are -not included.' However; - the present analysis 
can be used to examine the two-dimensional problem, and in particular, the 
sinusoidal pitching and plunging motion of an airfoil. The objective of this 
Appendix is to develop the equations for the unsteady pressure distribution 
associated with these motions. The resulting aerodynamic forces and moments 
can be' compiited from the zero and first order moments of the pressure 
distribution. 

The format for the analysis to be presented in this Appendix will 
parallel that of the main text. First a solution valid for small values of 
y (i.e., y = KM/(1-J^), where K ,= cob/u = Kb is the reduced frequency and O) is 
the circular frequency of .the .mptlpn) •will be derived by following the proce- 
dure developed in Appendix I. This will be followed by an analysis based on 
the theory in Appendix III which will be valid for large values of >. 


Solution for Small Positive Values of V 

The analysis presented in Appendix I showed that a solution for the 
noncirculatory flow field of order y is given in terms of a modified velocity 
potential which is a solution of Laplace's equation. This solution is 

given by Eq. (46) as 

= Z AnSinn-rye'"^ (108) 

n=i 


where the coefficients A^ are dependent on the blade motion (i.e., the A^'s 
are functions of the down'wash distribution). The value of these coefficients 
must be such that the fluid velocity normal to the surface of the airfoil is 
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equal to the velocity induced by the airfoil motion. The mathematical form 
of this bovindary condition is 


dXg “ at ax 


-b< X,<b, Xg = 0 (109) 


where h is the displacement of the airfoil from the X2 plane . The form of h 
for pitching motion about mid-chord is 


h = 


(no) 


while for plunging motion 


h 



( 111 ) 


where the parameters and h^ are the amplitudes of motion 

The value of the normal 
is given in terms of <|>^1) by 


The value of the noimal derivative of (b on the surface of the airfoil 
11) 


d<p _ yi -M= 
axg ’ b 


I 


? (I) 

-i/9M COS'*? -iwt a^ 


sin T7 


ac 


( 112 ) 


^^0 


whefe_ /3= k/ 1-1^. This resialt is derived from Eq. (l5)> by letting <x = 9 = 0 
and K U = <u. Thus from Eqs . (l08) and (ll2) 


_d±_ 

axo 




M' 


bsinT? 


^-ia,tg-i/3M"cos^ AnnsinT?! 


(113) 


which from Eqs. (109) and (llO) must be equal to (ioJO^^bX^^ + UaQ)e ^*^^for pitch' 
ing motion and from Eqs. (l09) and (ill) m\ist be equal to (-io)hQ)e“l‘^^ for 
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plunging motion. The coefficients Aj^ are evaluatefl hy employing the orthogo- 
nality relationship described in Eq. (5I) of Appendix I. The result for 
pltchl^ motion can he shown to he ' 


(a) 

An = 


Oob 


iwb 






20ob 


n-i 


(I - 




n > 


(nk) 


and for plunging motion the result is 



gicuhob i^~' 




n i I 


(115) 


The nonclrculatory pressure dlstrlhutlon on the surface of the airfoil 
is derived hy letting a b 0 ta 0 and K U ■ w in Eq. (2l) and comhlning the 
resulting expression With Eq. (108). The result is . 


s(i) 



Z A^ncosni? + i/3 Z Apsinn'^ 

' n»i ns| 


- Msecs’? (n 6 ) 


where the coefficients Aj^ are defined hy Eq. (ll4) or (ll5)* This nonclrcu- 
latory pressure dlstrlhutlon is unbounded at the airfoil leading and trailing 
edges . Only the trailing edge singularity must he removed . This is 
accomplished by the addition of a circulatory flow field. In Appendix I an 
analysis is presented which yields an expression for the circulatory pressure 
distribution on a swept airfoil encountering an oblique gust . The two-dimen- 
sional circulatory pressure distribution on the surface of the airfoil, 
obtained from the two-dimensional limit of that analysis, is 


s(2) _(2) -iySivi^cos'^ -iujt (117) 

p- = p e . . e 


51 



where 


_( 2 ) 

P 


^0^ C 
b sin^ 


[l +Zcos'^] 


(118) 


and where the parameter Z is defined hy Eq. (jl)- The constant C is assigned 
a value which will catise ps(2) 

at the trailing edge to be equal in magnitude 
but opposite in sign to p®(^) evalixated at the trailing edge. For plunging 
motion the value of C can be shown to be 


C = 


icuh^b 


(I - M^) 


[Jq{/3m 2) + iJ,(/3lV|2)]y/{l + Z) 


(119) 


while for pitching motion about mid-chord 


2 

Oob iw 


C = 


2{i-Mr 


( 120 ) 


I - 


[Jq{/3m2) + > / (I + Z) 


The pressure distribution on the surface of the airfoil for small values 
of y is obtained by combining Eqs . (ll6) and (llj)- The res\ilt is 



^i;A,ncosn77 + i^ EA^sinn^ 


-iwt 

e e 


C 

b sin'7 



( 121 ) 


The transfer functions associated with the unsteady motion of the airfoil 
are obtained by integrating the product of the pressure distribution and the 
appropriate weighting functions over the chord of the airfoil. For example, 
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the lift transfer function associated with plunging motion is dependent on 
the integral of the pressure distribution over the chord, and is given hy 
the equation 


T. 



U>y/l- 

ZttO^ 


io)t 

e 


( 222 ) 


where 




TT 


JS) 


sinT? d'’? 


2 0 


(123) 


and U2 = -iwh^. Substituting Eq. (l2l) into Eq. (l23) and performing the 
chordwise integration yields 





(22k) 


which reduces to 


C, = 


27T 1^ S ^ 2 .. >t 


L - , — n o 2 ^ nj. i/3Mn e 

-I " 


. 77 - (-iiuhj 




Jq(^ 8 M^) + iZJ,(/3M^) 


[Jq{/3M^) + iJ,(/SM^)] 


l + Z 


(125) 

“icjt 

e 


53 



when Eqs . (I15) and ( 119 ) are introduced. Thus the lift transfer function 
for plunging motion is equal to 



2i 


I - 


I nJp (/9 m2) 
n=l 


Jo(/ 3M2). + iZJ(/3M2)] r -I 

l + z ' I + iJ|(/3M2)J, , 


(126) 


The remaining transfer. functions can he derived hy a similar procedure. 

Solution for large Values of 7 

The analytical procedure presented in Appendix III can he used t o' develop 
the asymptotic form of the pressure distribution on a two-dimensional airfoil 
undergoing sinusoidal pitching or plunging motion. These expressions will 
approach the exact results as y become large . 


The asymptotic analysis will begin by considering the pressure field 
generated by a semi-indfinite plate oscillating sinusoidally in a pitching or 
plunging motion. The solution for the modified velocity potential associated 
with this pressure satisfies the two-dimensional, limit of Eq. (16) and the 
boundary condition 


dXg 



-i/ 3 iviS<i iwt 
e ' e. 


3h _3h_ 

at b aX| 


(127) 


on the surface of the plate. Equation-. (127) which is derived from Eqs. (15) 
and (109) is a .consequence of the boundary condition requiring the fluid 
velocity normal to the surface of the plate -to be equal to the velocity 
induced by the plate motion. The desired solution to the boundary value 
problem for can be obtained from an analysis based on the ¥einer-Hopf 
techn^ue, the essence of which can be found on pages i <-8 to 98 of Ref. 13 . 
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The pressure distrihution on the siirface of the plate is obtained by 
substituting the resulting solution for the modified potential into the 
two-dimensional limit of Eq. (2l). This procedure will yield 



for pitching motion about the leading edge (i.e., h = or^bX^ X^^O) and 






. KM 

( I +i ) e ' + M I 

'KM(I+M) y 27TX',“ 


for plunging motion (i.e., h = h^ As was shown in Appendix III the 

unsteady pressure distribution on a semi- inf inite plate approximates the 
pressure distribution in the neighborhood of the leading edge of an airfoil 
at high frequencies. These approximations (i.e., Eqs . (l28) and (l29))may ^ 
be extended to include the trailing edge region of an airfoil by combining 
the pressure field pk) with a field p'^) which is equal in magnitude but 
opposite in sign to p (1) downstream of the trailing edge (i.e., Xi 2 2, Xg = O) 
and has its normal gradient vanishing upstream of the trailing edge (i.e., 

Xg— 2, = O). The solution for such a pressure field can be obtained by 

following the procedure used in Appendix III. In that analysis an accelera- 
tion potential H is defined which is related to p'^' by the two-dimensional 
limit of Eq. (93) 




(130) 
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— f I f 

w]yere<Xj_ = Xj^ -2. Downstream of the trailing , edge , of the airfoil (i.e., 

Xi — 0, ^ =. O) the value of p'^'. is specified tp he eq^l.to the negative of 
p'^’. Thus H must’ equal 






(131) 


£n the boundary X^ — O, Xq = 0. The value of H on the upstream boundary 
X£-0, = 0 is given by Eq. (96) as 


( 132 ) 


H = 


TT 


- 3/2 ^ ( ■ >00 


/ M(U. xf;o) du X i 


X, < 0 , Xg'i 0 


where the function M(u, O) is defined by Eq. (97) and 


f()T;) = H 


■ X'Vo.'XgS 0’^' 


(133) 


For large values of y the function f may be approximated by the asymptotic 
form of H. This asymptotic expression is obtained by combining Eq. (l3l) 
with the asymptotic form of Eqs. (128) or (129). The result will be of the 
form 




5/ A 


J|V-, 


(I3i^) 


where 
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(l-i){l-M) 


^ 2(l+i).y2 + ^i^ I 


° ^ (KM)^/2 2yi + M y 27T (2 + Xf) My27rKM(l+M) ( 


- u^a , 
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My27rKM(|4M) 


2iKM 
I + M 


(136) 


and 


S,= 


KM' 


'■ I -m2 ’ 


8? = 


K M 
I -m2 


(137) 


for pitching motion about the leading edge. For plunging motion Eqs . (135) 
through ( 137 ) are replaced by 



U icuho 
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Substituting Eq. (13^) into Eq. (l33) and combining the result with Eq. (13&) 
yields 




||m(u,x/ ,( 


o)e"'^du X 


00 
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The evaluation of the integrals in Eq. (l^l) can again he performed hy the 
method of stationary phase (Ref. lU), which yields 


H - ( I - i) ^ g,(0) 


I + i 


- E 


y-(r+s,)x/ 


e ^ 92(0) 


+ i 




{/+ 82) x/ 


(li^2) 


X,' < 0 , X2 = 0 

Thus from Eq. (l30) the pressure distribution upstream of the trailing edge 
(i.e., X 2 = 0 ) is equal to 

,p( 2 > , p L.S.x/g |y-(r+8,)X,' ]| 

' (Ihz) 

Xi^+ wt] 

e 


where the functions and g 2 are defined by Eqs . (135) and ( 136 ) for pitching 
motion and by Eqs. ( 138 ) and (l39) for plunging motion. The parameters 5]_ and 
§2 are defined by Eqs. (13T) or (l^O). 

The addition of the pressure field p(2) to p(l) yields the asymptotic 
form of the pressure distribution p® on an oscillating two-dimensional airfoil 
in a compressible stream. The integral of the zero and the first order moment 
of this pressure distribution over the chord of the airfoil will yield the 
asymptotic forms of the unsteady lift and moment. 


+ e 


g^io) - E 


(r + S2)x, 



Limiting Behavior as K-+ 00 , M * 0 


The limiting form attained by the pressure distribution p® as the reduced 
frequency approaches infinity can be obtained by expanding the pressure 
distributions p^^) and p'^) in a power series of The result can be 
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shown to he dependent on the limiting form of p(l) , which is 

/oa[-iwaobx' + UQoj e'"^ 0< x|<2 (iM^) 


for pitching motion about the leading edge and 

ps=- p(i)2,^^a[-iwho]e-'"^ o<x/<2 ( 1 ^ 5 ) 


for plunging motion. Equations (144) and (l45) are precisely the results 
predicted hy piston theory (cf. Ref. 9) for pitching and plxmging motion. 

This result is to he expected, since piston theory may he used to predict 
the pressure distribution on an oscillating airfoil in a subsonic stream if 
the product of reduced frequency and ffech nvmher is large. The agreement 
between the limiting form of Eq. (lU2) or (1^3) and the results predicted 
by piston theory for pitching or plunging motion appears to be in conflict 
with the poor correlation between the two-dimensional limit of Eq. (36) and 
the result predicted by piston theory for the gust encounter problem. This 
apparent discrepancy can easily be resolved, by recalling that piston theory 
is only applicable to problems in which the local acceleration of the downwash 
distribution is much greater than the convected acceleration which is the 
case for the airfoil motion described above at high reduced frequency. For 
the gust encoimter problem the local acceleration and the convected accelera- 
tion of the downwash distribution are equal. 
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APPENDIX V 


NONPERIODIC TIME DEPENDENT FROBLEMS 


The objective of this analysis is to discuss the mathematical procedizre 
oy which the analysis presented in the main text can be extended to predict 
the unsteady airloads on an infinite swept airfoil encountering a nonperiodic 
gust. The downwash distribution resulting from the encounter of an infinite 
swept airfoil with a one -dimensional oblique gust is of the general form 


Ug = f [r - U cos ( 0 - a) t] 


Xp = 0 , -b<x<b 


(146) 


where 


r = Xjcosa + X 3 sina 


(147) 


This velocity field is convected at a speed of Ucos(^-a) in the r direction. 
The analysis in this Appendix assimies the existence of the Fourier transform 
in time of Eq. (l46). This transform is defined as 


I 


y 2 7T 'Leo 


I iwt , 

I Uge dt 


(M) 


while its inverse is defined to be 


I I 


.00, 


^2 = 


A I * io)t 


2 ir 


. I [) e aoj 


(149) 


The explicit form of the Fourier transform of Eq. (l46) is obtained by 
substituting Eq. (i4t) into Eq. (l48). The resulting integral may be evalu- 
ated by replacing the integration variable t with t* where r* = t - 


Ucos {d-otY 
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yielding 


ia;r 

"I _ ucos(0-a) 


f(oj) 


-b<x,<b, X2=o 


(150) 


where 


f(w) = 



(151) 


The spatial dependence of Eq. (150) is identical to that for a one-dimensional 
sinusoidal gust (i.e.-, Eq. (lO)). The remaining equations to he transformed 
are Eq.- (t) for the; pressure field and Eq. (8) for velocity po^e^tial. These 
eqToations are transformed hy multiplying each term hy e and inte- 

grating the result with respect to time from minus to ;^us infinity. This 
procedure leads to the result 


P = -P^ 


+ UCOS0 


dX, 


+ Lisin© 


dX, 


(152) 


for the Fourier 


transform of Eq. 


(T) and 


o'" I - M ^cos^© - M‘'sin‘©-^ 


d 


.2_. 2, 


d> 


dX 


dx; 


- 2a‘ 


_.M^ ^ + M^cosesine^ 

a dx, a aXa dx.dx. 


(153) 


cu 


<P = 0 


for the Foxirier transform of Eq. ( 8 ). Equation (l 53 ) can he reduced to the 
Hemholtz equation (i.e., Eq. (16)) for a modified velocity potential This 
transformation is accomplished hy introducing dimensionless coordinates 
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(i.e., Eqs . (ll), (l2), and (13)) and the following form for <f> 


« A 


i 

Xg) e 


-/3M^COS0X| + 


Kb sinaXj 
(I - 


into Eq. (l53) which yields 


axf 


a^<l> 

axf 


? A 

y^<J> = 


0 


(I5i^) 


(155) 


The variables yS and T which appear in Eqs. (l5^) and (l55) are defined by 
Eqs. (iH) and (l^) of the main text while the variable K is related to w by 
the equation 


K 


O) 

Ucos(0- a) 


(156) 


The boundary conditions for can be explicitly expressed in terms of 4 and 
its spatial derivatives . The transformation of the downwash boundary 
condition, Eq. (9)^ is 


A 

axg 



X 2 = 0, -b < X| < b 


(15T) 


which, when combined with Eqs. (1^9) and (l52), yields 


a<i» 

axg 


- M^cos^e 


f(w)e 


i/3x, 


-I < X,< I, Xg = 0 (^58) 


The transformation of the Kutta-Joukowski condition requires p® to vanish 
downstream of the trailing edge of the airfoil. The mathematical form of 
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this condition is obtained by combining Eqs . (152 ) and (15^)> yielding 


f\ s 
p 



UcosQ 

b 


aX| 


i /3<t 


-i/9M®COS^0X| + 

e 


iKbsinaXj 
(I -M^cos^©)''^ 


0 (159) 


X, = I, X 2 = 0 


The field equation (l55) and the boundary conditions in Eqs. (I 58 ) and (l59) 
form a boundary value problem for 4 > which is identical to that for the 
sinusoidal gust discussed in the main body of text. Thus, the solution for 
the modified velocity potential associated with the sinusoidal gust problem 
can be used in conjunction with Eqs. (l52) and (l5^) to predict the Fourier 
transform of the press\are distribution p® on the surface of the airfoil. 

The integrals of the zero and first order moments of the Fourier transform 
of the pressure coefficient induced by the gust are equal to 


and 


Cl = 


2 7Tf C0S8 

I - M^COS^ e 




KsinaXj 


( 160 ) 


A 

Cm = 


i RsinaXj 


7Tfcos9 { 'Ksinax^i 

22 ‘ 


/i-M cos 0 


( 161 ) 


which are the Fourier transforms of the lift and moment coefficients. The 
functional forms for the variables T^ and Tj^^ are defined in the main text. 
Finally the lift and moment coefficients resulting from the encounter of the 
airfoil with the gust are obtained by inverting the Fourier transforms, 

Eqs. ( 160 ) and ( 161 ), yielding the results 

(162) 


C|_ XjSina - Ucos(6- a)t] = 


cos0cos(9- a) iK[x 3 sina -ucos(0-a)r 
7 - . a 

- M^COS^S 
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and 


( 163 ) 


C|^[x 3 sina - Ucos(0- a)t] = 


y cos0cos(0-a) 




^ \ - M^COS^0 '00 




iK XjSina 


■ UCOS(0-a)t 


dK 


A typical application of the analysis in this Appendix is to the prohlem 
of predicting the unsteady air loads generated hy a helicopter rotor blade 
encountering a tip vortex from a previous blade . This problem can be analyzed 
by considering the interaction between an infinite swept airfoil and a two- 
dimensional rectalinear vortex which has an induced dovmwash field given by 

- r{r-Ucos(^- g)t| , 

^2 ■ / 2 2 ) 

27t < [r - Ucos(0-a)t] + h ^ 

where h is the vertical position of the vortex beneath the X2 plane and 
r is the strength of the vortex. The Fourier transform of Eq. (l 64 ) is 


,1 

U2 


sgn(K).y/^ 


iKr 

ir/2Tre 


ucos (0-a) 


-|R|h 

e 


(165) 


Thus from Eq. (15O) 


f(K) = sgn(K) 



-IRIh 
i r/27re 

u cos (0- a) 


(166) 


Introducing Eq. (166) into Eqs . (162) and (163) yields 


C. = 




rcos 0 i fCO _ iR [xjsina - Ucos( 0 - a)t] _ (167) 

• sgn(K)e T. e dK 


M^COS ^0 
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for the lift coefficient and 




r COS© i 


8U^ 


r/i2 2q 

M COS 6 


rCO _ "l^lh 

-/ sgn(K)e T e 

-00 


iKiXjSina - ucos(fl-a)t] 


dK 


(168) 


for the moment coefficient. Numerical results for the lift or moment coef- 
ficient can he obtained hy introducing the appropriate transfer function into 
these equations and integrating the resulting expression numerically. This 
procedure has been used to obtain the results presented in the main text. 
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APPENDIX VI 


CQMPtPTER PLOW DIAGRAM FOR COMPUTING THE AIRLOADS GENERATED 
BY THE ENCOUNTER OF A HELICOPTER ROTOR BLADE WITH A TIP VORTEX 


The flow diagram for the computer program that was developed to predict 
the aerodynamic loads generated hy the encounter of a helicopter rotor "blade 
with a tip vortex is presented in this Appendix, along with a list of input 
and output variables . 


Input Varia"bles 

1. Azimuth position at which the encounter occiirs . 

2. Advance ratio. 

3- Rotational tip Mach number, 

if. Angle of encounter. 

5 • Rotor blade semi-chord . 

6. Vertical displacement of the vortex beneath the plane of the rotor. 
7- Radial position at which the encounter occurs. 

Output Variables 

1. Mach number at the position at which the encounter occurs. 

2 . Sweep angle . 

3- Time history of the nomialized lift coefficient. 

h. Time history of the normalized moment coefficient. 
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Evaluation of Eqs . 

(167) and (168) by 
Numerical Integration 


Output 
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FIGURE 1. COORDINATE SYSTEM GEOMETRY 
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EDGE 


FIGURE 2. GEOMETRY OF THE ENCOUNTER 
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FIGURE 3. GEOMETRY DEFINING THE RELATIVE MACH NUMBER M* 
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EXACT (THEORY OF REF. 7) 
• APPROXIMATE 



FIGURE 4. A COMPARISON BETWEEN EXACT AND' APPROXIMATE 
EXPRESSIONS FOR THE LIFT TRANSFER FUNCTION 
FOR 7 > AND MACH NUMBER 0.3 






EXACT (THEORY OF REF. 7) 



FIGURE 5. A COMPARISON BETWEEN EXACT AND APPROXIMATE 
EXPRESSIONS FOR THE LIFT TRANSFER FUNICTION 

FOR 7 > 0 AND MACH NUMBER 0.6 ; 
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EXACT (THEORY OF REF. 7) 
APPROXIMATE 



FIGURE 6 ._ _A COMPARISON BETWEEN EXACT AND APPROXIMATE 
EXPRESSrONS FOR THE MQM^J TRANS FER FUNCTION 
FOR 7>6 AND MACH NUMBER 0.3 
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EXACT (THEORY OF REF. 7) 



FIGURE?. A COMPARISON BETWEEN EXACT AND APPROXIMATE 
EXPRESSIONS FOR THE MOMENT TRANSFER FUNCTION 
FOR y >0 AND MACH NUMBER 0.6 




APPROXIMATE 



FIGURE 8. A COMPARISON BETWEEN EXACT AND APPROXIMATE 
EXPRESSIONS FOR THE LIFT TRANSFER jFUNCTION 
FOR 7 = 0 AND MACH NUMBER 0.3 





EXACT 

APPROXIMATE 



FIGURE 9. A COMPARISON BETWEEN EXACT AND APPROXIMATE 
EXPRESSIONS FOR fHE LIFT TRANSFER FUNCTION 
FOR 7 = 0 AND MACH NUMBER 0.6 









REALTivi,^^ 


FIGURE 11. A COMPARISON BETWEEN EXACT AND APPROXIMATE 


EXPRESSIONS, FOR THE MOMENT TRANSFER FUNCTION 
FOR 7 = 0 AND MACH NUMBER 0.6 




REAL Tl 


note I„, Tl = 0 
EXACT (THEORY OF REF. 7) 
APPROXIMATE (EQ. (34) ) 
APPROXIMATE (EQ. 39) ) 



GUSt WAVE NUMBER,' Kb 


FIGURE 12. A COMPARISON BETWEEN EXACT AND APPROXIMATE 
EXPRESSION FOR THE LIFT TRANSFER FUNCTION 
FOR 7^ < 0 AND MACH NUMBER 0.3 
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REAL T|_ 



FIGURE 13. A COMPARISON BETWEEN EXACT AND APPROXIMATE 
EXPRESSION FOR THE LIFT TRANSFER FUNCTION 
FOR 7^ < 0 AND MACH NUMBER 0.6 
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REAL 


note = 0 

EXACT (THEORY OF REF. 7) 
APPROXIMATE (EQ. (35) ) 
APPROXIMATE (EQ. (40) ) 


0.3 X 10~^ 



figure 14. A COMPARISON BETWEEN EXACT AND APPROXIMATE 
EXPRESSIONS FOR THE MOMENT TRANSFER FUNCTION 
; FOR 7 ^ < 0 AND MACH NUMBER 0;3 
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FIQU RE 15. A COMPARISON BETWEEN EXACT AND APPROXIMATE i 
EXPRESSIONS FOR THE MOMENT TRANSFER FUNCTION 
FOR 7 ^ < 0 AND MACH NUMBER 0.6 




LIFTING LINE 



FIGURE 16. LIFT COEFFICIENT FOR A PARALLEL ENCOUNTER 



MOMENT ABOUT OUARTER CHORD 



FIGURE 17. MOMENT COEFFICIENT FOR A PARALLEL ENCOUNTER 





CONFIGURATION AT 



FIGURE 18. LIFT COEFFICIENT FOR AN OBLIOUE ENCOUNTER 





APPROACH MACH NO. 0.6 
MOMENT ABOUT QUARTER CHORD 



FIGURE 19. MOMENT COEFFICIENT FOR AN OBLIQUE ENCOUNTER 


GEOMETRIC CONFIGURATION 
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FIGURE 20. LIFT COEFFICIENT PER UNIT SPAN AS A FUNCTION OF TIME 
AND VORTEX VERTICAL DISPLACEMENT FOR A PARALLEL ENCOUNTER 





APPROACH MACH NO. = 0.6 



FIGURE 21. MOMENT COEFFICIENT PER UNIT SPAN ABOUT THE QUARTER-CHORD 
AS A FUNCTION OF TIME AND ' VORTEX VERTICAL DISPLACEMENT FOR A PARALLEL ENCOUNTER 





CONFIGURATION AT TIME t 
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NORMALIZED LIFT COEFFICIENT 


Rj - ROTOR RADIUS 


0.85 Ry 

0.753 Rj 

0.674 R j 


CONFIGURATION AT TIME t = 0 



FIGURE 24. NORMALIZED LIFT DISTRIBUTION FOR SECOND-QUADRANT VORTEX 

ENCOUNTER OF A UH-ID ROTOR AT 90 KT 
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NORMALIZED MOMENT COEFFICIENT, 


Rj - ROTOR RADIUS 


0.85 Rj 
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0.674 Rj 



X3 

NONDIMENSIONLESS SPANWISE COORDINATE, Z = — g— sin a 


FIGURE 25. 'WORMALIZED MOMENT DISTRIBUTION FOR SECOND-QUADRANT VORTEX 

ENCOUNTER lOF' A UH-ID ROTOR AT 90 KT 
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